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THE HITTING TIMES WITH TABOO 
FOR A RANDOM WALK ON AN INTEGER LATTICE 



Ekaterina VI. Bulinskayauo 



Abstract 

For a symmetric, homogeneous and irreducible random walk on Z d , d E N, having 
zero mean and a finite variance of jumps, we study the passage times (taking values in 
[0,oo]) determined by the starting point x, the hitting state y and the taboo state z. 
We find the probability that these passages times are finite and analyze the tails of their 
cumulative distribution functions. In particular, it turns out that for the random walk 
on Z d , except for a simple (nearest neighbor) random walk on Z, the order of the tail 
decrease is specified by dimension d only. In contrast, for a simple random walk on Z, 
the asymptotic properties of hitting times with taboo essentially depend on the mutual 
Q^ location of the points x, y and z. These problems originated in our recent study of 

Ph ' branching random walk on Z with a single source of branching. 

^-> . 
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O ■ 1 Introduction 

A random walk is a classical model being the source of numerous interesting problems having 
elegant solutions. The monographs [IT], [19] and [22] devoted to various properties of random 
walks have become the reference books for many researchers. The recent publications such 
as [12] and [18] demonstrate its non-vanishing popularity. The reason is the abundance of 
applications of this simple probabilistic model (see, e.g., [3] and [28]). Moreover, a lot of new 
complicated models were constructed on the basis of random walk (one can mention those 
investigated in jTT] and [20]). Thus, many assertions of the theory of random walks have been 
obtained as auxiliary (technical) lemmas needed for understanding other models of interest. 

The present paper also provides results concerning random walk on integer lattice Z d (d G N) 
which arose during the study of another model, namely, branching random walk (BRW) on Z d 
with a single source of branching (different modifications of the model were considered, e.g., in 
papers [2], [9], [25] and [30]). Its main features are the following. If the starting point of the 
process is x G Z d , x 7^ 0, then the parent particle performs random walk on Z d until hitting the 
origin. If the starting point is 0, or the parent particle has just hit it, then after exponentially 
distributed time the particle may die producing a random number of offsprings. Otherwise, 
it leaves the origin (jumps randomly to a point x' G Z d \{0}) and behaves afterwards as a 
memoryless particle starting at x'. At the birth moment the newborn particles are located 
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at the origin. They and all their descendants evolve according to the scheme described above 
independently of each other as well as of their parents history. 

In the framework of BRW a natural question arises: what is the limit distribution, as 
t —7- oo, of the number of particles /i(£; y) at each point y G Z d at time £? Evidently, fi(t; y) 
essentially depends on the trajectory of the parent particle. Indeed, if the parent particle hits 
the origin earlier than the time t then //(£; y) can take any value in Z+. Otherwise, //(£; y) — 1 
or fi(t; y) — if the parent particle is located at time t at point y or outside it, respectively. 

This is the motivation to introduce the notion of hitting time with taboo for a random walk 
on Z d . More precisely, for any x,y,z G Z d such that y ^ z, let function H Xt y tZ (t), t > 0, be 
(improper) cumulative distribution function (c.d.f.) of the first hitting time (or the first return 
time if x — y) of point y if the starting point of the random walk is x and the point z is a 
taboo state. In the next section we give a formal definition of this function. In the present 
work we find the limit value H x ^^ z (oo) = lim^oo H Xtyz (t) and analyze the asymptotic behavior 
of H Xj y Z (oo) — H xyjZ (t), as £ — > oo, for a symmetric, homogeneous, irreducible random walk on 
7L d [d G N) having zero mean and a finite variance of jumps. 

It turns out that the most interesting case is d = 1 since there are two quite different 
kinds of asymptotic behavior of H Xty>z (t), t — > oo, depending on whether the random walk is 
simple (nearest neighbor) or not. It is also worth mentioning that, for the random walk on Z d , 
except for a simple random walk on Z, the value H Xy y jZ (oo) G (0, 1) and the order of decrease 
of function H Xjyz (oo) — H x<yz (t) is determined by dimension d only, regardless of x, y and z. In 
contrast, for a simple random walk on Z, the mutual location of points x, y and z determines 
the value H Xy y tZ (po) G [0, 1] as well as the order of decrease of H x ^ jZ (oo) — H Xt y iZ (t) as t — )■ oo. 

Finally, we recall that the properties of hitting times (or, more generally, passage times 
including first entrance and last exit times) and taboo probabilities for a Markov chain with 
stationary transition probabilities were exposed in [T3] and [23], respectively. For the taboo 
probabilities see also [31] and references therein. The counterpart of the function H x>yjZ (t) for 
a Markov chain with stationary transition probabilities has been used earlier (see, e.g., [10J, 



p. 202, and [16J, p. 31). However, it was treated as an auxiliary tool and has not been studied 
per se. 

2 Main results 

We assume that all random variables (taking values in [—00,00]) are defined on a probability 
space (Q, J 7 , P). We study a class of random walks on Z d more general than that proposed in 
118], Ch.l, Sec.2. 



Definition 1 For d G N, let A = (a(x,y)) xyeZ d be a matrix such that elements a(x,y) > 
whenever x 7^ y, a(x,x) < 0, X]«ez d a ( x ' v) = ® an( ^ su Pxez d |a(:r,:r)| < 00. Then a continuous 
time Markov chain S = {S(t),t > 0} having the state space Z d , i.e. S(t) : fl — > Z d for each 
t > 0, and generator A is called a random walk on Z d . 



The existence of such a Markov chain is guarantied by Theorem 2 in [T5], Ch.3, Sec.2. 

From here on we consider only a symmetric, homogeneous, irreducible random walk S 
having a finite variance of jumps. Symmetry and homogeneity mean that for x, y G Z d one 
has a(x, y) = a(y, x) and a(x, y) = a(0, y — x) =: a(y — x), respectively. Random walk is called 



irreducible if for any x,y G Z d there exists t > such that P(5*(t) = y\S(0) = x) > 0, i.e. all 
the points of the lattice Z d can be reached. Furthermore, a random walk is said to have a finite 
variance of jumps if XLez d ll 3 -!! 2 * 3 ^ #) < oo where is the origin of Z d and || ■ || is a norm in 
this space. 

The conditions imposed on elements a(x, y) allow us to use an explicit construction of the 
random walk on 7L d with generator A (see, e.g., Theorem 1.2 in [5], Ch.9, Seel). According 
to this construction S is a regular jump process with right continuous trajectories. A jump of 
the process from state x G Z d to state y G Z d , x ^ y, occurs with probability a(x,y)/a(x,x) 
independently of the earlier evolution of 5*. Moreover, for transition times of the process 
r(°) := and r^ := inf{£ > r^ 1 ) : S(t) ^ ^(r^- 1 ))}, n > 1, the following statement 
holds. Random variables {r( n+1 ) — t^}^ =0 are independent and each of them has exponential 
distribution with parameter a := — a(0,0). In what follows we consider the version of the 
process 5* constructed in this way. 

Note that for the sequence {t^"''}'^ =1 of transition times of 5* one can consider embedded 
chain {S n ,n G Z + } defined by way of S n := S(t^). Mention in passing that {S n ,n G Z + } is 
a discrete-time homogeneous Markov chain having transition probabilities a(x,y)/a for x ^ y, 
x,y G Z d (see, e.g., [5], Ch.8, Sec.4). 

Recall some other facts concerning the introduced random walk on 7L d . Let p(t; x, y), t > 0, 
i,t/6 Z d , be the transition probabilities of the random walk, that is, 

p(t; x, y) := P(S(t + u) = y\S(u) = x) = P(S(t) = y\S(0) = x) 

for any u > 0. Then for h — >■ 0+ 

p(h;x,y) = a(x,y)h + o(h), if x^|/, 
p(h;x,x) = 1 + a(x,x)h + o{h). 

Due to [5], Ch.8, Sec. 3, the transition probabilities of the process S* satisfy the backward 
Kolmogorov equations 

d~T~ = *l2 z&d a ( x > z ")P( t > z >y^ p(0;x,y) = 5 y (x). (1) 

As usual, S y (x) is equal to 1 or for x = y or x ^ y, respectively. The Fourier transform 
applied to (JTJ) yields (see [25] and [29], Ch.2, Seel) 

p(t;x,y) = -±- d [ e mt+i(0,y-,) ddj xyeZ * ( 2 ) 

where (•, ■) stands for the scalar product in Euclidean space R d and 

= V" o(0, z)cos(z,6), 9 G [-7i,' 



7r| 



Furthermore, by the method of steepest descent one derives from ([2]) the asymptotic relation 

v(t;x,y)~—, x,yGZ, t^oo. (3) 



Here 

1 ( d' 



(27r)"/V|det^(0)| ™ V ' V W. 
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A similar result for transition probabilities of a discrete time random walk on Z d was established 
in [22], Ch.2, Sec.7. 

Moreover, according to [29J, Ch.2, Seel, equality (J2]) implies that 

p(t; 0, 0) - p{t; x, y) ~ ™ "^ , x,yeZ d , t -> oo, (4) 

where 

7 d (z) : = _L^ /" (^) 2 e §W*(°)^) d v , z g Z rf . 



2(2tt)^, ; 

Set G\(x,y) := / °° e p(t;x,y) dt, A > 0, x,y G Z d , that is, G\(x,y) is the Laplace 
transform of the transition probability p(-; x, y). In view of (J2J) the introduced function can be 
represented in the form 

1 r f i{S,y-x) 

Ga(x,1/) = — -^/ t— r^de. (5) 

(27T) d ./r-ff-id A - 



Taking into account ([3]) it is easily seen that for d > 3 the Green's function Go{x,y) is 
finite for all x,y G Z°\ However, for d = 1 or d — 2, one has lim.\^o+ G\(x,y) = oo. This 
phenomenon corresponds to transience of the random walk for d > 3 and its recurrence for 
d — 1 or d — 2 (for more details see, e.g., [18], Ch.4, Seel and 2). 

By virtue of (TjJ the function \im\^ 0+ (G \(0 , 0) — G\(x,y)) is finite for all d G N and 
x, y G Z d . Therefore we can define 

/% f a lim(Gi(0,0)-Gi(0,x)), if X / 0, 

Ml):= \i if , = 0. 

Equality (J5J) allows us to give another formula for Pd(x) if x 7^ 

a f cos(x,9) — 1 ,. 

Now we are able to introduce some basic notation and formulate our main results. Set 
t := inf{£ > : S(t) 7^ S(0)}, i.e. the stopping time r (with respect to the natural filtration of 
the process S) is the time of the first exit from the starting point of the random walk. Taking 
into account the explicit construction of the random walk one has r = r^ almost surely (a.s.) 
and, consequently, G(t) : = P(r < t\S(0) = x) = 1 - e~ a \ t > 0. 

Let T y>z '■= inf{£ > r : S(t) = y,S(u) 7^ z, r < u < t] for y, z G Z d , y 7^ z. As usual, 
inf{£ G 0} = 00. We call the stopping time r y>z a hitting time of the state y with the taboo state 
z. Denote its (improper) c.d.f. given the starting point x by H x>yjZ (t) := P{T y>z < t\S(0) = x), 
x,y,z G Z d , y 7^ z, t > 0. 

Note that the random variable r~ := inf{£ > : S(t + r) = y, S(u) 7^ z,t < u < t + t}, 
y,z £ Z , y 7^ z, is not a stopping time w.r.t. the natural filtration of the process S. However, 
T y;Z can be also called a hitting time of the state y with the taboo state z, since in this case we 
count time not from the moment of the process start but from the moment of the first exit out 



of the starting point. In a similar way, we set H xyz {t) := P{r yz < t\S{fS) = x), x,y,z G Z d , 
y^z,t>0. 

Obviously, r y>z = T~ + r a.s. Moreover, in view of strong Markov property of the random 
walk w.r.t. the stopping time r (see, e.g., [5], Ch.8, Sec. 4) the random variables t~ z and r are 
independent. Hence, T y ^ z has an absolutely continuous (improper) distribution as a sum of two 
independent random variables, one of them having a density. On the other hand, the explicit 
construction of S shows that the distribution of r~ z has an atom at zero when x ^ y, namely, 

H x, y A°) = a{x,y)/a. 

For x, y G Z d such that y 7^ 0, we specify the functions Cd(x, y) by formulae 

_ pxiy - x) + pi{x) - pi(y) airy 2 'fi(pi(y - x) - pi(x) - pi(y)) 2a-nxy^\ 
4a7T7! pj(y) p^y) 

n ( x p2Jy - x) + p 2 (x) -p 2 {y) 

C2{x,y) := , 

4a7 2 

^a(pd(y -x)+ p d (x) - p d (y)) 

a(d-2){G (0,0) + G {0,y)f 

The following theorems are the main results of the paper. 

Theorem 1 Let x,y,z G Z d be such that y ^ z. Then for the random walk S on 7L d , except 
for a simple random walk on Z, one has 

tt ( \ Pd( x ~z)+ Pd{y -z)- p d (y -x) 

H X! y iZ (oo) = — -, r G (0, 1), d = l or d = 2, (7) 

2p d (y - z) 

tt r^\ g o(°» °)pd(y ~ z )~ Gq(0, 0)p d (y -x) + G {y, z)p d (x - z) 

H x ,y, z {oo) = tt G 0, 1 , d > 3, (8) 

Pd{y - z){G {0, 0) + G (y, z)) 

Moreover, as t — > oo, 

iW°°) - H x>y , z (t) ~ Cifr-^-*) for d = lf (9) 



^ / \ _ g ld\Hd\y _j>j-r ^d w _ Hd^yj) , . Q 



C 2 (x - z,y- 


-z) 


ln£ 
C d (x-z,y- 


-z) 



H X:y>z (oo) - H x>y>z (t) v ln ^ '- for d = 2, (10) 

H x>y>z (oo) - H x , y , z (t) ~ ^"^if"^ /^ rf>3 (11) 

where C d (-,-), d G N, are positive functions defined above. 

Theorem 2 Let S be a simple random walk on Z. // x,y,z G Z and y ^ z then 

v 2 7/ — X 

1 - H X} y !Z (t) ~ — r , x <y < z or z <y < x, (12) 



'07T 



v/t 



^12/ — ^| \J2aix\/t 

x — z 

— H Xty>z (t) = o(e" a£< ), 2<x<y or ?/ < x < z, (14) 

2/ ^ 

-^-iJ z ,^(t) = o(e- a£< ), x = z, (15) 

H x , y , z (t) =0, x<z<y or y < z < x, (16) 



/or some e G (0, 1). 



Theorem 3 Theorems [T] and [2] hold true if function H XyV)Z at the left-hand sides of fl7|)- f|T6l) is 
replaced by H~ , whereas the right-hand sides of these formulae remain intact. 

The proofs of Theorems [T] and [2] are given in Sections 4 and 5, respectively. They bear on 
application of the Laplace-Stieltjes and the Laplace transforms to the main equations involving 
H x ,y,z and H x ,z,y as well as functions H xy , H zy , H xz and H yz . Here H XtV (t), x,y G Z d , t > 0, 
is the (improper) c.d.f. of the hitting time of point y given the starting point x. After that we 
employ the Tauberian theorems techniques and, for d sufficiently large, we essentially use Faa 
di Bruno's formula for the n-th derivative of two functions superposition, n > 1. An essential 
difficulty arising on this way is to prove the positivity of functions Cd(x,y), x,y G Z d , j/^ 0, 
under the conditions of Theorem [TJ Theorem [3] proved in Section 6 can be viewed as a corollary 
of Theorems [1] and EJ 

3 Auxiliary results 

For a nonnegative function n(t) and a nonnegative nondecreasing function x{t)-> t > 0, we set 

/•OO /"OO 

k(\):= e- xt K(t)dt, x(\):= e - xt d X (t), A > 0, 

Jo Jo 

whenever the integrals exist. Recall a useful relation linking the Laplace-Stieltjes transform of 
the bounded function x(t) and the Laplace transform of x(oo)— xif) where x(°°) := h m t->oo x(j), 
namely, 

(X(oo)-XJ(A) = - , A > 0. (17) 

This relation results from the formula of integration by parts. 

Introduce a stopping time T y := inf {t > r : S(t) = y} called the hitting time of point y G Z d . 
Denote its (conditional) c.d.f. given the starting point x by H xy (t) := P(r y < t\S(0) = x), 
i,l/6 Z d , t > 0. The limit value H XtV (oo) := linn^oo H Xty (t) and the asymptotic behavior of 
function H xy (oo) — H x>y (t), as t — > <x>, were found in papers [H] [E], [25] and [27] for y = and 
x G Z d . In Lemma [1] we reformulate these results in a more general way. 



Lemma 1 Let x,y G Z d . Then 



H X;y (oo) = 1 /or d — 1 and d = 2, (1£ 

1 
aG (0,0) 

Go^y) 

G (0,0) 



-f^xx(oo) = l — — for d > 3 and V = x, (19) 

-G (0,0) ■/ - y v ; 

g»,»(°°) = ^7^ f° r d >3 and V^ x - (2°) 



Moreover, as t — >■ oo, one /ias 



l-ff^(t)~ Pl(y ^ /or d = l, (21) 



a7i7r 



v^ 



l-iJ^(t)^ P2(?/ f /or d = 2, (22) 

a 72 Int 



PROOF. In view of definitions of r and r x as well as of conditional probability and a Markov 
chain we see that 

p{t- x, x) = P(S(t) =x,t> t\S(0) = x) + P(S(t) = x,t< t\S(0) = x) 
= P(r > t\S{0) = x) + P{S(t) =x,t x < t\S{0) = x) 

= 1 - G(t) + I P(S{t) = x\t x = u, S(0) = x) dP(r x < u\S(0) = x) 
Jo 

= 1 - G(t) + [ P(S(t) = x\S(u) = x) dH XtX (u). 
Jo 

Thus we deduce the following formula 

p(t;x,x) = 1-G(t)+ / p(t-u;x,x)dH x>x (u), xeZ d , t > 0. (24) 

Jo 

Note that (JSJ implies the identity G\(x,x) = G\(0, 0) for any x G Z d . Therefore, application 
of the Laplace-Stieltjes transform to ( 124)) leads to 

^- (A) = * ~ G A (z,x) = X " G A (0,0) = ^ ' o(A) - (25) 

In particular, by the Laplace-Stieltjes transform uniqueness (see, e.g., [H], Ch.13, Seel) it 
follows that 

H x>x {t) = H , {t), t > 0, x e 7L d . (26) 

For different d, asymptotic properties of H 0t0 (t), t — > oo, were established in [6], [25] and 
[27] . More exactly, the limit value Hq (oo) := \im. t ^oo Hq (t) and the asymptotic behavior 
of H 00 (oo) — H 00 (t), t — > oo, were found for the function Hq (J:), t > 0, related to //o,o by 
Ho,o(t) = G* H 00 (t) (as usual, * denotes the convolution). However, these results are easily 
extended to the case of function H 00 (t). Hence, due to (|2"6"|) relations ffT8|) . (fT9l and f[2Tj) (T23|) 
are proved for x = y. 

Let x j^ y. Combined definitions of r y , conditional probability and a Markov chain imply 



p(t;x,y) = P(S(t) = y,r y <t\S(0) = x) 

= f P(S(t) = y\r y = u, S(0) = x) dP(r y < u\S(0) = x) 
Jo 

= f P(S(t)=y\S(u) = y)dH Xiy (u). 
Jo 

Consequently, for x ^ y and t > we get the formula 

p(t;x,y)= / p(t-u;y,y)dH x , y (u). (27) 

Jo 

By virtue of (jHJ) identities G\(x, y) = G\(x — y, 0) = G\(y — x, 0) and G\(y, y) = G A (0, 0) are 
valid. Then application of the Laplace-Stieltjes transform to (1271) yields 

ir m G x {x,y) G x (x-y,0) G x (y-x,0) - - 

H ^ X) = G^y) ~ G A (0,0) " G A (0,0) " H ^ o{X) ~ F ^ o(A) " (28) 

7 



Now using the Laplace-Stieltjes transform uniqueness we obtain 

H x>y (t) = H x _ y>0 (t) = H y _ x>0 (t), t > 0. (29) 

The limit value H z ^ (oo), z G Z d \{0}, and the asymptotic behavior of H z ^(oo) — H Zj0 (t), 
£ ^ oo, were found in [7] and [8]. It follows that by ((29]) relations flTg), (f20]) and (l2T])-(l2"3"j) are 
established for x ^ y. The proof of Lemma [T] is completed. 

The next lemma provides a linear integral equation involving functions H x ^ yz (t) and H x ^ z ^(t) 
which is the basis of our theorems proofs. 

Lemma 2 Let x,y, z G Z d and y ^ z. Then the following equation holds 

H x , y {t) = H x ^ z {t) + [ H z , y (t - u) dH x ^ y (u), t > 0. (30) 

Jo 

Proof. Firstly note that the random variables r y and t VjZ coincide a.s. on event {r VtZ < oo} 
in view of their definition, that is T y \[r y ^ z < oo) = T y ^ z \{r y ^ z < oo) a.s. (here 1(B) stands for the 
indicator of a set B). Consequently, {r y < t,T y<z < t} = {t VjZ < t] for each t > 0. Moreover, 
one can see that r y l(r Z)y < oo) > T z>y I(T Z!y < oo) a.s. Therefore, for every nonnegative t, one 
has {r y < t,T yjZ > t} = {r y < t,r yjZ = oo} = {r y < t,T z>y < oo} = {r y < t,r z ^ y < t}. These 
relations allow us to write 

H x , y (t) = P(r y < t, r y>z < t\S(0) =x) + P(r y < t, T y , z > t\S(0) = x) 
= P(r y , z < t\S(0) =x) + P(r y < t, r z , y < t\S(0) = x) 



H x , y , z (t) + / P(r y < t\T Zyy = u, 5(0) = x) dH x>z>y (u). 
Jo 



Using definitions of r y and r z<y as well as that of a Markov chain it can be verified that 
P( r j/ < A T z,y — u , S(0) = x) = P{j y < t — u\S(0) = z). The previous reasoning supplemented 
with this equality entails the desired statement. Lemma [2] is proved. 

The obtained equation has a natural interpretation. Namely, the term at the left-hand side 
of ( 13 Op accounts for the random walk trajectories starting at point x and hitting point y until 
time t. Every path of such kind belongs to one of two types. The trajectories of the first type 
do not pass point z before hitting y. They are taken into account by the first summand at 
the right-hand side of ( 130]) . As to a trajectory of the second type, it hits point z at time u, 
< u < t, before reaching y so that the part of the trajectory after hitting z is a path starting 
at z and reaching y until time t — u. The second summand at the right-hand side of (13"0"|) is 
responsible for the trajectories of the second type. 

The next statement is a Lemma |2] corollary which provides the Laplace-Stieltjes transform 
of a solution of equation (|30|) . 

Corollary 1 If x,y,z G Z d and y ^ z then for any A > 

H x , y , z (X) = ^ (A) 7^ (A) ^ (A) . (31) 

1 - H Zt y(\)Hy, z (\) 



Proof. Due to Lemma [2] we have a system of two linear integral equations in the functions 
H x>y>z (t) and H x>z>y (t) 

±lx,y\b) = Hx,y,z\t) ~r tl x ,z,y * ^z,y\p)i 

H X)Z {t) = H XtZt y(t) + Hx^^* Hy tZ (t). 

Application of the Laplace-Stieltjes transform to each equation of the system leads to a new 
system of algebraic equations in H x ,y,z(ty and H XtZ>y (X) 

H X .y(X) = H Xt y iZ [\) + H XiZi y[X)H Zt y[\), 

H x . z {\) = H XtZ>y (X) + H Xt y !Z (X)Hy >z (X). 

Solving this system we obtain (I3T]) . The proof is complete. 

The next proposition clarifies the invariance of functions H Xyy>z {-) indexed by x,y,z G Z d , 
z j^ y, under shifts and reflections of space Z d . 

Corollary 2 Let r, x,y,z G Z d and y 7^ z. Then 

H-x+r,y-\-r,z+r\t) ^x,y,z\^J ■> tl — x,— y,— z\t) ^x,y,zV') j " — "■ v / 

Proof. The desired assertion is a direct consequence of formulae ( 1261) . ( 129]) and ( 13T1) . 

For proving Theorem 1 we also need the following lemma. Its derivation is based on con- 
sidering of all the possible jumps of the random walk after latency period r. 

Lemma 3 For x,y,z G Z d such that y ^ z, one has 

H- ytZ (oo) - H- ytZ (t) = J2 ^fl(H r , y , z (oo)-H r , y , z (t)), t>0. (33) 

r£Z, d ,r^x,rj^y,rj^z 

Proof. In view of definitions of random variables r and r~ one has 

y> z 

H- y , z (oo) - H- yz (t) = P(t<r- z <oo\S(0) = x) 

E P(*<r^<oo,5(T) = r|5(0) = ar) 

rd1i d ,r(f:{x,y,z} 

J2 P(S{T) = r\S(0) = x)P(t<T- z <oc\S{T) = r,S(0)=x). 

r£Z d ,r$:{x,y,z} 

Here P(£ < r~ z < oo\S(t) = r, 5(0) = x) = P(t < r y>2 < oo\S(0) = r) by the identity 
T y,z = t + Ty^ z a.s. and due to strong Markov property of the random walk w.r.t. the stopping 
time r (see, e.g., [5], Ch.8, Sec. 4). Moreover, according to the explicit construction of S the 
relation P(S(r) = r\S(0) = x) = a(x,r)/a is also true. Combining the obtained equalities with 
the previous argument one comes to (I331) . Lemma [3] is proved. 

The next result which is known as Faa di Bruno's formula plays an important role while 
establishing Theorem [T] for d > 3. 

Theorem 4 (see [24J, Appendix, Subsec.9) Let functions W (u) andV(u) have n-th deriva- 
tives. Then for the n-th derivative of function U(u) = W(V(u)) the following formula is valid 

"*•'<»> = £ w( —^^rr. ©"' ® " (£)" ' ' - < 34 > 

na,n2,n3,...eZ_|-: x ' ' 

ni+2n 2 +3n 3 + ...=n 



The last auxiliary result will be useful for proving Theorem [2j 
Lemma 4 For each x G N the following equality holds true 



1 — cos 
Proof. For each n G N one has 



' ■■™- rB dB = 2*x. (35) 



cos n9 — cos (n + 1)9 f cos (n — 1)9 — cos n9 

d9 — / d9 



1 — cos 9 / 1 — cos 9 



2cosn9 — cos (n + 1)9 — cos (n — 1)9 f 2cosn9(l — cos9) . 

d9 = d9 = / 2cosn9d9 = 0. 



1 — cos 6* 7 1 — cos 

— 7T 

Whence for all n G N the values of the following integrals are the same and equal, for instance, 
the value of the integral for n — 1 

* cos (n- 1)0- cos n0 f/1 /"* 1 - cos ,„ 

at) = - rf0 = 27T. 



1 — cos ./_*■ 1 — COS 

Hence, for each lENwe obtain 

/*" 1 — COS.T0 v-^ f n COS (n — 1)0 — COS7Z0 v^ 
— = > / d9 = > Z7T = Z7TX. 
^ 1-COS0 ^7_ 1-COS0 £^ 

Thus, equality (135]) is established and the proof of Lemma @] is complete. 

4 Proof of Theorem H 

Firstly, in view of Corollary [2] it suffices to prove Theorem [T] for z = 0. Secondly, we have to 
consider separately the cases d = 1, d = 2 and d > 3, since the demonstration of the theorem 
essentially depends on dimension d. Before passing to the non-simple random walk on Z let us 
write down some formulae which have the same form for all deN. 

Let x,y G Z d \{0} and x ^ y. Using (fTTJ) , flU, (J2SD and flS} as well as recalling that 
J °° e _A4 (l — G(t)) (it = (A + a) -1 , A > 0, we deduce the following relations 

(H (Z?T n \(\\ H **s>(°°) G x (x,y)G x (0,0) -G x (0,x)G x (0,y) 
(H x , y , (oc) - H x , y , )(X) = — (G2(0,0)-G*(0,y))A ' (36) 

W»)-W(A) - ~ x { Gl(0,0)-Gl(0,y))(\ + a)\> (3?) 

W«»)-WW " ~ x - X + { Gl(0,0)-Gl(0,y))(\ + a)\- (38) 
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4.1 The case d = 1 

In this subsection we focus on a non-simple random walk on Z although some arguments are 
valid for a simple random walk on Z d as well. At the beginning we find the limit value H x>yj o(oo) 
for x, y G Z, y ^ 0, due to the formula H Xj y$(oc) = limA_+o+ H x ,y,o(ty ■ To this end we write 
the following asymptotic representation of the function if ,r(A), r G Z, 

g 0r (A) = l- Pl(r) ^? + o(v / A), A^0+, (39) 

071V 71 " 



which results, in view of (I17p . from application of Tauberian theorem (Theorem 4 in [T4], Ch.13, 
Sec. 5) to relation (121"]) . Then substituting ( 139]) in (13T]) and taking into account (129]) we calculate 
limA_s>o+ H Xt y }0 (\) and, consequently, we get formula for H Xj y t0 (oo) which is true for a simple 
random walk on Z as well. However, to complete the derivation of (J7J) for d = 1 we have to 
show that H X) y j0 (oo) G (0, 1). Note that by virtue of ( 118]) equality ( 130]) implies the equivalences 
H x , yt o(oo) = 43- H X)0t y(oo) = 1 and H Xjy>0 (oo) = 1 <^> H X)0>y (oo) = 0. Hence, to verify that 
H x ,y t o(oo) G (0, 1) for all x, y G Z, y ^ 0, it suffices to check that H Xt y j0 (oo) > for all x, y G Z, 
y 7^ 0. In its turn the latter holds true if, for instance, relation (J9]) is satisfied with Ci(x, y) > 
for all i,|/6Z,!/^0. Thus, the initial problem of validating that H Xt y t0 (oo) G (0, 1) is reduced 
to proving ([9]) and showing positivity of the function Ci(-, •). The rest part of the subsection 
is devoted to establishing the last two claims. 

The demonstration of is based on formulae (]3"oT) . (I3"7j) and (13"%]) in which we substitute 
asymptotic decompositions of functions G\(0,0) and G\(0,0) — G\(0,r), as A — > 0+, for 
r = x, r = y and r = y — x. Due to Tauberian theorem (Theorem 4 in [14J, Ch.13, Sec. 5) the 
asymptotic behavior of G\(0, 0) for d = 1 is encoded in relation ([3]), namely, 



G A (0,0) = ^ + o(-y , A^0 + . (40) 

Furthermore, using the same Tauberian theorem we infer that 

G A (0,0)-G A (0,r)=a- 1 pi(^)-2v / ^7iWv / A + o(v / A), A -+ 0+, r G Z\{0}. (41) 
Indeed, a" 1 pi(r) = J °° (p(t; 0, 0) — p(t] 0, r)) d£ < oo on account of (jl]) and 

/•oo 

G A (0,0)-G A (0,r) aT l Pl (r) = e- xt (p(t; 0, 0) - p(t; 0,r)) d£ - a"Vi(r) 

JO 

/OO / ft \ />oo 

e_A * ( / (p(u;0,0)-p(u;0,r))du\ dt- X e~ Xt a~ l p x {r) dt 

/OO / /"00 \ 

e" A M / (p(«;0,0)-p(u;0,r))d«J dt. (42) 

Here J°° (p(w; 0, 0) — p(u; 0, r)) dw ~ 271 (r) t -1 / 2 (as £ — > 00) by formula (jl]) and Theorem 31 in 
[26J. Now we substitute relation (J4T]) and the formula for H x>yt0 (oo) appearing in (J7J) into (136]) . 
( 137]) and ( 138]) . After collecting terms we substitute ( 140]) in the obtained formulae. Omitting 
tiresome calculations we deduce that for x, y G Z, y 7^ 0, and A — > 0+ 

(«,,,(^ff„„.)(A) ~ ^ I '+^-'>-"-M 

4 a v^Tiv A 
Qy 7 ^ (pi(y)7i(x) + p x (y - xffxjy) - pi{x)^i(y) - pi{y)ii(y - x)) 

pt(y)V\ 
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This can be rewritten in the form 

(M?^(^^|^. A^0+, (43) 

given that 71 (r) = 7rr 2/ jf, r G Z. The latter identity is valid by virtue of definitions of 71 and 



7i(r) since 71 = 1/y — 27T0"(O) and 



r 2 /"+ 00 „ r 2 /•+ 00 , r 2 



!-./-.' ' " -4tt0 // (O) v /-0 // (O) 7-oo" " -20" (0) ^-271-0" (0)' 



—00 



We took into account that the variables change u = v^/ — 0"(O) reduces the integration to 
writing J_ °° u 2 e~ u I 2 du = y/2n. Note that relation (J4"3l is valid for a simple random walk on 
Z as well. However, in this case the function Ci(x,y) takes nonnegative values (not strictly 
positive values). 

If Ci(x,y) > then application of Tauberian theorem (Theorem 4 in [14], Ch.13, Sec. 5) to 
relation (143j) leads to the desired statement ([9]). Thus, to complete the proof of Theorem 1 for 
d — 1 we only need to verify the positivity of function Ci(x, y), x, y G Z, y 7^ 0. 

According to Corollary [2] we may assume that y > 0. Then there exist three possible relative 
positions of points x, y and 0, namely, x > y, < x < y and x < 0. Let us consider at first 
the case x > y. To check strict positivity of Ci(x,y) for such x and y we search for a lower 
estimate of it. To this end note that pi(x) + p±(y — x) — pi(y) > in view of formula for 
Hy tX fi(oo) appearing in (J7j) and the evident inequality H ytXt0 (oo) < 1. Obviously, C\(0,y) > 
and, consequently, p\{y) — Aa 2 T: 2 y 2 ^\ > 0. Combining these inequalities we see that for x > y 

n( v (pi(x) + Pi(y - x) - pi{y)){pj{y) - Aa 2 n 2 y 2 ^) 2 Pl (y - x)a^y 2 ^l 

4a7rp 2 (y)7 1 p{{y) 

2an^ly{x - y) > 2p x {y - x)a-ny 2 i[ 2am*f\y{x - y) 

pi(y) ~ Pi(v) pi(y) 

Before passing to the case < x < y we derive one more useful relation. Recall that 
T y ,o > r y0 a.s. because t Vi0 = r y0 + r a.s. Since the random variable r is exponentially 
distributed, P{r y0 < 00) = P(t^ < 00). Therefore P(t < T y> o < 00) > P(t < r yQ < 00) and 
hence H Xt y t0 (oo) — H x>yt0 (t) > H~ y0 (oo) — H~ (t). Due to Lemma and the monotonicity 
property of the Laplace transform it follows that 

(H x>yi0 (^ H Xjy , )(\) > J2 ^^-(Hr^o&h 1 H r ,y,o)(X), A > 0. (44) 

Now turn to the case < x < y. We show the positivity of Ci(x,y) successively for 
x = y — 1, x = y — 2, . . ., x = 0. li x = y — 1 then there exists a point r > y such that 
a(y — 1, r) > as otherwise the random walk is simple. Therefore, by virtue of (14*31) and of the 
just established inequality Ci(r, y) > 0, < y < r, there is a summand of the order 1/a/A at 
the right-hand side of ( 1441 . Consequently, by ( 1431) the left-hand side of ( )44l) is also of the order 
1/vA, A — t- 0+, and the relation Ci(y — 1, y) > with y > is proved. Recall that the random 
walk under consideration is irreducible, that is, the greatest common divisor of all r G Z\{0} 
such that a(r) > equals 1 (see [4], Ch.12, Sec. 3). This implies for x = y — 2 that there exists 
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a point r>y or r = y — 1 such that a(y — 2, r) > 0. Since for such r we have just shown 
the positivity of Ci(r,y), in view of (j43p it means that the left-hand side of ( l4"4"j) has the order 
of growth 1/a/A, as A —)• 0. So we conclude that C\(y — 2,y) > with y > 0. Considering 
successively the cases x = y — 3, x = y — 4, . . ., x — 0, we prove according to the described 
scheme that Cx(x, y) > for all i£ [0,j/),iGZ. 

The case x < < y is also treated by indirect methods, since the immediate estimation of 
C\(x, y) for such x and y is rather difficult. We consider in succession x = — 1, a; = —2, .... At 
first we assume that for the random walk S there exist at least three points r\ > 0, r 2 > and 
r3 > such that a(ri) > 0, a{r-i) > and 0(7*3) > 0. If cc = — 1 then at the right-hand side of 
( T44j) there is at least one summand indexed by r where r > 0, r ^ y and o(— 1, r) > 0. Similarly 
to the previous discussion we infer that the left-hand side of ( l44"j) has the order of growth 1/vA, 
as A — > 0+, and Ci(—l,y) > 0. If x = — 2 then at the right-hand side of (j4"4"|) there exists at 
least one summand indexed by r = — 1 or r > 0, such that r /?/ and o(— 2, r) > 0. Since for 
such r the positivity of Ci(r,y) has been already proved, in view of (l4"3"j) and (j4"4"l) we deduce 
that Ci(— 2, y) > 0. The positivity of Ci(x,y) for a* = —3, x = —4, ... is verified in the same 
manner. 

We continue to deal with the case x < < y assuming that there exist exactly two points 
ri > and r 2 > such that a(r*i) > and 0(7*2) > 0. Let x = — 1. If there is a point r > 0, 
r ^ y, such that a(— 1, r) > (i.e. t*i = r + 1 or r 2 = r + 1) then by virtue of (1431) and (j4"4"|) as well 
as of the verified inequality Ci(r,y) > we get the desired relation C\{— l,y) > 0. Otherwise, 
one has r\ — 1 and r 2 = y + 1. If y > 1 then a(— 2, j/ — 1) = a(y + 1) > and, consequently, 
C\{— 2,y) > in view of (j4"3"|) . (jHJ) and the established earlier estimate Ci(y — 1, y) > 0. 
Therefore, if r\ = 1, r 2 = y + 1 and ?/> 1 we obtain Ci(— 1, y) > due to (j4"3~|) and (jHJ) as well 
as the relation a(— 1, —2) = o(l) > and the just proved inequality Ci(— 2, y) > 0. However, 
for ri = 1, r 2 = y + 1 and y = 1, the previous arguments do not work and we have to employ 
the following facts concerning the embedded chain {S n ,n e Z + }. 

For the embedded chain {S n ,n G Z + } introduce a passage time T := min{n > : S n > — 1}. 
According to [I], Theorems 15.1 and 15.2, 

P(T < cx)|5(0) = -1) = 1, P(T > 7i|5(0) = -1) ~ 4=, n^oo, (45) 



where c is some positive constant. 

For brevity we write P-i(-) := P(-|5(0) = —1). Since we assume that r± — 1 and r 2 = 2, the 
process {S n , n G Z + } can perform jumps to the nearest-neighbor points or to the next nearest 
points. Hence, 

P_i(T = n, S n = -1) = P_x(T = 7i, S^i = -2, ^ n = -1) + P_i(T = n, S„_! = -3, S n = -1) 

= P_!(5 fc < -1, < k < n - 1, S n _x = -2, 5 n = -1) 
+ P^(S k < -1, < k < n - 1, 5 n _! = -3, S n = -1) 

= a( ~ 2 '~ 1) P-i(g fc < -1, < A; < n - 1, S n _! = -2) 
a 

+ ^~ 3 '~^ P-i(g fc < -1, < k < n - 1, Sn-i = -3) 

> min{a(1) ' a(2)} P_ 1 (r = n). (46) 
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Let us return to a lower estimate of i?_i t i,o(cx>) — H-i^ t o(t) = P-i(t < ri 7 o < oo), as t — > oo. 
Denote by iV = {N(t),t > 0} the Poisson process constructed by means of the random sequence 
| r (™+ 1 ) _ r^l^Q, i.e. iV is the Poisson process with intensity a. Considering all the possible 
jumps of 5 1 , taking into account (|46p and invoking the explicit construction of S we derive that 

oo 

P-i(t<Ti, <oo) = 5^P_i(ri, = r (n) ,t < n, < oo) 

n=l 
oo 

= 53 P_x(5 fc ^ 0, 3 k ± 1, < A; < n, S n = 1, t < r (n) < oo) 

n=l 
oo 

= 53 P-iOSfe < -1, < k < n - 1, 5 n _i = -1, 5 n = 1, JV(t) < n) 

n=l 
oo 

= J3 P_!(5 fc < -1, < A; < n - 1, S n _i = -1, 5 n = l)P(JV(t) < n) 



n=l 



> Q( 1,1; ^ p -i( 5 fc < -1, < Jfe < n - 1, 5 n _i =: -l)P(iV(t) < n) 

n=l 



■oo 



" (2) J] P_!(T = n - 1, SU =. -l)P(iV(t) < n) 



a 

n=l 



n-1 



> AJ n ;, i ;J 53P-i(T = n-i)53 



-att„4-\k 



a(2)min{a(l),a(2)}^ n .^e" ai (at) 



a 2 ^-f "^ '^ fc! 

n=l fe=0 

a(2)min{a(l),a(2)} y, p i(T>fe) g^W (4?) 

fc=0 

Thus, finding a lower estimate of P_i(t < 7~i : o < oo) is reduced to establishing a lower estimate 

of the last series. It is not difficult to check that (N(t) — at)/\/at — > £ with £ ~ A/"(0, 1), as 
t — y oo. Consequently, P(N(t) G (at — by/at, at + b\/al)) — >■ P(£ G (—6,6)) > 0, as t — )• oo, for 
each 6 > 0. Hence, in view of (1451) 



oo —at( +\k 

V P_i(T > fc)- — -^- > P(T > [at + 6Vat]) P(-W(£) G (at - 6\/at, at + 6^)) 

fc=0 

P(£e (-&,&)), t^oo. (4? 



at 

Combining ( I4U1) . ( T4T1) and ( 1481) as well as taking into account the monotonicity property of the 
Laplace transform we conclude that C\{— 1, 1) > 0. 

Thus, we have shown that C\(— l,y) > for each y > 0. Let us pass to the starting point 
x = —2. Recall that we assume that there exist exactly two points r\ > and r 2 > such that 
a(ri) > and a(r2) > 0. If there is point r > 0, r ^ y, or r = — 1 such that a(— 2, r) > (that is, 
r\ = r + 2 or r 2 = r + 2) then due to ( 1431) and (14"4"|) as well as the verified inequality C\(r, y) > 
we get the desired estimate C\{— 2,y) > 0. Otherwise, one has r± — 2 and r 2 = 2/ + 2. In 
particular, a(— 4, —2) > and a(— 4, y — 2) > 0. Since the random walk is irreducible, one can 
guarantee that y — 2 ^ 0. Consequently, Ci(— 4, y) > in view of ( 1431) . ( j44l) and the earlier 
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established inequality C\{y — 2, y) > with y > and y ^ 2. In its turn C\{— 2, y) > by virtue 
of (J43]) . (jUJ and the just proved estimate C\{— 4, y) > 0. Therefore, the positivity of Ci(x, y) 
for x = —2 and y > is demonstrated. For x = —3, x = —4, . . ., verification of Ci(x,y) 
positivity, y > 0, is implemented similarly to the case x = —2. Consequently, Theorem [T] for 
d — 1 is proved. 

4.2 The case d = 2 

As in the previous subsection one can find the limit value H x ^ y $(oc) for x,y G Z 2 , y 7^ 0, by 
the formula H Xt y j0 (oo) = liniA-s.o+ H Xt y ;0 (X). For this purpose we write the following asymptotic 
representation of the function H 0r (\), r G Z 2 , 

# 0>r (A) = 1 + -&£L + o (J-V A^0+, (49) 

a 72 In A \ In A / 

which results from application of Tauberian theorem (Theorem 4 in [14], Ch.13, Sec. 5) to 
relation (T221 and from employing (fT7|) . Substituting (j4"9l in (J3"T|) and taking into account (T2"9l 
we find lim A _ s . + H x ^ y0 (X) and, thus, H Xt y t0 (oo) is obtained. However, to complete the derivation 
of (j7|) for d = 2 one has to show that H Xty0 (oo) G (0, 1). Similarly to Subsection 4.1, equalities 
( TTBl) and ( 1301) entail that H Xjyj o(oo) = <^> H Xt o ty (oo) = 1 and H Xtyi o(oo) = 1 <^> H Xj o, y {oo) = 0. 
Therefore, for establishing that H x>yt o(oo) G (0, 1) for all x, y G Z 2 , y 7^ 0, it suffices to verify 
the inequality H X) y t0 (oo) > for all x, y G Z 2 , y ^ 0. In its turn the latter holds true if, 
for instance, relation fTTOT) is satisfied with C 2 (x,y) > for all x,y G Z 2 , ?/ ^ 0. Hence, the 
initial problem of checking that H x<yt0 (oo) G (0, 1) is reduced to proving of (TTU1) and showing 
the positivity of function C 2 (-, •). The rest part of this subsection is devoted to validating the 
last two claims. 

The proof of ffTUj) bears on formulae (J36|) - fl38|) in which we substitute asymptotic decom- 
positions of functions Ga(0,0) and Ga(0,0) — G\(0,r), as A — > 0+, for r = x, r = y and 
r = y — x. Due to Tauberian theorem (Theorem 2 in [TJ], Ch.13, Sec. 5) the asymptotic 
behavior of G\(0, 0) for d = 2 is implied by ([3]), namely, 

G A (0,0) = -7 2 lnA + o(lnA), A^0 + . (50) 

Using the same Tauberian theorem and relation (J4"2]) . valid for all d G N, we deduce that 

G A (0,0)-G A (0,r) = a-V2(r)+72(r)AlnA + o(AlnA), X -)• 0+, r G Z 2 \{0}. (51) 



Indeed, in view of fll]) one has a _1 p2( r ) — J"^ (p(£; 0, 0) — p(t; 0, r)) dt < 00 and by formula (J3J 
and Theorem 31 in [26J one infers that f. (p(u; 0, 0) — p(u; 0, r)) c?m ~ 7 2 (r)/£, as £ — >• 00. Now 
we substitute relation floTj) and the formula for H Xt y y0 (oo) appearing in ([7]) into (136]) - (138]) . After 
collecting terms we substitute (1501 in the obtained formulae. Omitting tedious calculations we 
establish that 

(H Xi j, i0 (oo) - f/a^oXA) — — , A -> 0+, z,j/ G Z , yf=0. 

—4 o 72 A In A 

If 6*2(0;, y) > then application of Tauberian theorem (Theorem 4 in [13], Ch.13, Sec. 5) to the 
last asymptotic relation gives ( FlOl . Thus, to complete the proof of Theorem [1] for d = 2we only 
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have to verify the positivity of C 2 (x, y) for all x, y G Z 2 , y 7^ 0. To prove the claim we employ 
the method of ad absurdum. Assume that there exist 2 and y such that x G Z 2 , y G Z 2 \{0} 
and 6*2(2, y) = 0. Then the definition of 6*2(2, y) implies that 

P2(x) + p 2 (y ~ x) - p 2 (y) = 0. (52) 

In view of (J7J) the last equality is equivalent to the relation H yjXj0 (oc) = 1. Therefore, formulae 
(TT8|) and (I30I) entail that H yt o tX (t) = and H yyXy0 (t) = H VtX {t) for £ > 0. Hence, by virtue of 
(|22|) one has G2(y, 2) > and 6*2(2/, x ) — p 2 (y — x) / (a j 2 ) . Taking into account the definition of 
C2(y, 2) the last equality can be rewritten as p 2 (y) — p 2 (x) = 3p 2 (y — x). This relation combined 
with assumption ( 1521) leads to contradiction. Hence, C 2 (x,y) > for all x,y G Z 2 , y 7^ 0. The 
proof of Theorem [T] for d = 2 is complete. 

4.3 The case d > 3 

In this subsection we assume that d > 3. Basing on the formula H x ^ )0 (oo) = lim A _> + H xy $(\), 
2, y G Z d , y 7^ 0, we find the limit value H x ^ )0 (oc) with the help of CHI . (J2UJ) and (13"T1) . To prove 
(JH) completely we have to check that H x ^ (oo) G (0, 1). Obviously, in view of (fT9l) . fl20l) and 
( |30|) the relation H Xj y t0 (oo) < H Xty (oo) < 1 holds true. For verifying inequality H x ^ (oo) > it 
suffices to establish ffTTj) and show the positivity of C d (x, y). So the rest part of the subsection 
is devoted to validating these two assertions. 

Further on we suppose that 2 and y are fixed points within Z d \{0}. Formulae (I3"6"j) - (1551) 
can be rewritten in the following form 

(u n^T~ u \(w pd(y)(pd(y - 2) + p d (x) - Pd (y)) 

(H x , yt0 (oo) - H x , y , )(\) = a 2 (Go(0; o ) + Go(0;?/)) -M* °> V) + -M* *> V) + J 3(A; x, y) 

G {0,0){p d {y-x)+p d {x)- Pd {y)) a(G (0, 0)G (2, y) - G {0, x)G {0, y)) 

-J4a;2/J — - J 5 {A;y), 



a(G (0,0) + G (0,y)) ^ '* y Pd(y)(G (0, 0) + G (0, y)) 

(53) 

(M»)-W(A) = a2(Go(0 , 0) + G (0, y)) Jl(A; " y) - a(G (0, 0) + G (0, y)) ^ y) 

(;,,(0 "' /) -J 5 (A;y) + ij 6 (A;y,y), (54) 



p d (y)(G (0,0) + G (0,y)) °^ a/ a 

(^o(^P^o)(A) = aa(Go(0> 5 + ^(0^)) Jl(A; °' y) - a (Go( G 0°) ( + J o( 0,y)) J4(A;?/) 

+ r vr ^TXr rn ^ (A;y) - -J 6 (A;0,y) (55) 

Pd(y)(Go(0,0) + Go(0,y)) a 

where for r G Z d and A > 

^ d (A;r) (tf d (A;y-r)-# d (A;y))G A (0,0) 

Ji(A;r,y) := J 2 {X;r,y) •- 



G 2 (0,0)-G 2 (0,y)' -v>'^- G 2 (0,0)-G 2 (0,y) 

m x _ (K d (\;0)-K d (\;r))G x (0,y) _ K d (X; 0) - K d (X;y) 

31 ' "' " J •" G 2 (0, 0) - G 2 (0, y) ' MA > V) - G\ (0, 0) - G 2 (0, y) ' 

K d (A;0)-^(A;y) G A (0,r) 

J s(A; y) := „ , n n , — „ , n v , ^e(A; r, y) :-- 



G A (0,0)-G A (0,y)' ^ ' '^- (A + a)(G 2 (0,0)-G 2 (0,y)) 
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and 

G (O,r)-G A (O,r 
K d (X;r) := 



A 

Let us briefly describe the further scheme of the proof of ffTTl) . At first we differentiate 
[d/2] — 1 times (symbol [•] stands for the integer part of a number) the both sides of (153"]) . 
( 15"4"|) and ( l5T)j) w.r.t. variable A at point A > 0. Then we retrieve the asymptotic behavior 
of function j{ (A; r, y), as A — > 0+, for each r G Z d . Moreover, we verify that the func- 

tio^4 [d ^\X;r,y), jf ^(Ajr.v), df ^(A;*,), J^H^y) and jf/ 2 ^(A; r, y) are 
°(^i (^> r ' ?/))> as A — )• 0+, for each r. Afterwards one can employ Corollary 43 in J2H] and, 

thus, find the asymptotic behavior of the functions H x> y t0 (oo) — H xy0 (t), H 0t y >0 (oo) — H Qy0 {t) 
and H y ^fi(oc) — H y} y^ (t), as t — > oo. 

Before studying the limiting behavior of j{ (A; r, y) we have to investigate the asymp- 

totic properties of the functions K d (X;r), (G A (0,0) - G x (0,y)y\ (G|(0, 0) - G 2 x (0,y))-\ 
G A (0, r)/(G|(0, 0) - Gl(0, y)), G A (0, r)/(A + a) and their derivatives up to the order [d/2] - 1 
for each r G Z d . 

Integrating by parts it is easy to check the following equality 



oo 

-At 



K d (X;r)= e (/ p(u;0,r)du\ dt, r G Z d , A > 0. (56) 

Then for each n G N function K d (X; r) has the n-th derivative w.r.t. A > which can be 
calculated by the formula 



/oo / />oo \ 

e' M t n I / p(u; 0, r) du J dt 



(57) 



for each r G Z d . For convenience we set K d (A;r) := K d (X;r). Let us analyze for each r the 
asymptotic behavior of the n-th derivative (0 < n < [d/2] — 1, n G Z) of the function K d (X; r) 
at point A, as A —7- 0+. To this end we need the relation implied by (|3]) and Theorem 31 in 
|26j . namely, 

/OO Q 

p(«;0,r)d«~ (d _ 2 w/ 2 -i » *->°°> rGZ "- ( 58 ) 

Then in view of (1571) one has, as A — >• 0+, 



/OO / POD \ 

t n f / p(w;0,r)dw] d£<oo, 0<n< 



d 



1. (59) 



Moreover, using Tauberian theorem (Theorem 2 in [T4], Ch.13, Sec. 5) we observe that, as 
A^0+, 

(_l)rf/2-2^/2-2) (A . r) _ 27^ ^1 onlyforeven d (60) 

(ci — 2) A 

Note that for uneven d the numbers [(d — l)/2] — 1 and [d/2] — 1 coincide whereas for even 
d one has [(d — l)/2] — 1 = d/2 — 2. Employing relations ( |5"7|) and ( 158]) as well as Tauberian 
theorem (Theorem 4 in [T3], Ch.13, Sec. 5) we infer that, as A — > 0+, 

( _ 1)W2| - 1/ , r ,- 1 . (A;r) f^. « - —en, 

l P^A» lf rf 1SeVen - 
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We stress that the right-hand sides of (160)) and (|6T|) do not depend on point r G Z d . 

Now we investigate the asymptotic properties of function G A (0, 0) — G A (0, y) and its deriva- 
tives up to the order [d/2] — 1, as A — ¥ 0+. Recall that for A > one has the identity 
G A (0,0) — G\(0,y) = / °° e~ xt (p(t; 0, 0) — p(t;0,y))dt where the asymptotic behavior of the 
function p(t; 0, 0) — p(t; 0, y) (as t — > oo) is given by (0J. Then by virtue of the counterpart of 
( 157)) for the function G A (0, 0) — G A (0, y) we get that, as A — > 0+, 



/■oo 

4) n (G A (0,0)-G A (0,j/))M-> / t n (p(t;0,0)-p(t;0,y))dt<oo, 0<n< 

Jo 



-1. (62) 



Using Theorem H] we are ready to estimate the asymptotic growth of the n-th derivative, 
< n < [d/2] - 1, of the function (G\(Q,Q) - Ga(0,2/)) _1 , as A ->■ 0+. For this purpose we 
consider function V~ x as the external function W(V) and we substitute G A (0,0) — G\(0,y) 
instead of internal function V(-) in TheoremHl Since W^(V) = iy~ l y' = (— l) n n\ V~ n ~ l for 
every n e Z + and G A (0, 0) — G A (0, y) — > a^ 1 p d (y) > 0, as A — >• 0+, then in view of Theorem 0] 
and relation (16"2~)) we obtain the estimate we are interested in 



G x (0,0)-G x (0,y) 



(«) 



0(1), A^0+, 0<n< 



- 1. 



(63) 



Let us turn to the study of the asymptotic behavior of the function G A (0,r) and its n-th 
derivatives, < n < [d/2] — 1, as A — > 0+, for each r G Z d . In view of the counterpart of ( 1571) 
for G A (0, r) as well as of relation (J3J) we see that, as A — > 0+, 



r(n) , 



-l) n G A nj (0,r)^ / t n p(t;0,r)dt < oo, 0<n< 



d-1 



- 1. 



(64) 



Moreover, employing Tauberian theorem (Theorem 2 in [2], Ch.13, Sec. 5) we find that 



-l) d/2 - 1 Gf 2 - 1) (0,r)~ 7d ln^ A 



0+, only for even d. 



(65) 



As mentioned above, the numbers [(d — l)/2] — 1 and [d/2] — 1 coincide for uneven d whereas 
for even d one has [(d - l)/2] - 1 = d/2 - 2. 

Now for estimating the asymptotic behavior of the n-th order derivatives, < n < [d/2] — 1, 
of the function G A (0, r)/(A + a), as A — > 0+, it suffices to apply the Leibniz formula as well as 
relations (1541 and (1651). Indeed, 



GA(0,r; 
A + a 



(n) 



J2c k n G{ k \0, 



k=0 



A + a 



(n-fe) 



Y,C k n G{ k \0, 



k=0 



(-l) n - k (n-k)\ 
(A + a) 



n— fc+1 



where n G Z + and, hence, for each r G Z d we obtain the desired estimates 



G x (0,r) 



(n) 



A + a 

GxC0 1 r)^ (d/2 - 1) 

A + a 



0(1), A^0+, 0<n< 



d- 1 



-1, 



O In 



A 



A — t- 0+, only for even d. 



(66) 
(67) 
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Repeating the argument used for deriving relations ( I64p and ( 165]) we find the asymptotic 
behavior of the function G A (0,0) + G A (0,y) and its derivatives up to the order [d/2] — 1, as 
A^0+, 



d-1 



-l) n (G A (0, 0) + G A (0, y)) (n) -» J £>(£; 0, 0) + p(t; 0, y))dt < oc, < n < 

o 

(_l)^-i (GA(0) 0) + G A (0, 2 /)) W2 - 1) ~27 d lni only for even d. 



1, (68) 



(69) 



To investigate the asymptotic properties of the function (G A (0,0) + G A (0,y)) _1 and its 
derivatives of order n, < n < [d/2] — 1, as A — > 0+, we employ Theorem H] once again. 
We take 7 _1 as the external function 1^(7) and consider the function G A (0,0) + G A (0,y) 
as V(-) in Theorem H Since W^(V) = (y- 1 ) 1 ^ = (-l) n n! 7~ n - 1 for all n G Z+ and 
G A (0, 0) + G A (0, ?/) ->■ G (0, 0) + G (0, y) > 0, as A ->• 0+, due to Theorem 1 as well as re- 
lations (168]) and ( 169]) we get the estimates 



1 



G A (0,0) + G A (0,y) 



(n) 



0(1), A^0+, < n < 



d- 1 



(70) 



(d/2-1) 



O ( In - ) , A -)- 0+, only for even d. (71) 



G x (0,0)+G x (0,y) 

To establish asymptotic behavior of the function (G\(0, 0) — G x (0, y)) _1 and its derivatives 
up to the order [d/2] — 1, as A — > 0+, it suffices to apply the Leibniz formula along with 
estimates (155]) . (17D 1 and (171 )1 . Indeed, 



(n) 



% G^(0,0)- 05(0,1/), 

for n G Z+, and, consequently, 

1 



£<« 



w 



G A (0, 0) - G A (0, y) ; VC A (0, 0) + G A (0, y) 



1 



(n-fc) 



G*(0,0)-G*(0,y) 



(n) 



0(1), A^0+, < n < 



d-1 



-1, 



(d/2-1) 



O In 



A — > 0+, only for even d. 



(72) 



(73) 



G|(0,0)-G|(0,y)y ~ V"A / 

Let us investigate the asymptotic properties of the function G A (0, r)/(G|(0, 0) — G\(0, y)) 
and its n-th derivatives, < n < [d/2] — 1, for each r G Z d , as A — > 0+. The Leibniz formula 
implies that for all A > and n G Z + 



G x (0,r) 



G|(0,0)-Gl(0,y) 



(n) 



Ew(°- 



fc=0 



G*(0,0)-G*(0,y) 



(n-fc) 



Taking into account ( 164"]) . ( 165]) . ( 1721 and ( 1731 . for each r G Z d , we infer that 



G x (0,r) 



Gl(0,0)-Gl(0,y) 



(n) 



0(1), A^0+, 0<n< 
19 



d-1 



(74) 



G\(0,r) \ (d/24) / l 

O I In — ) , A — > 0+, only for even d. (75) 



Gl(0,0)-Gl(0,y)J V A, 

Since all the basic relations are established, we turn to finding the asymptotic behavior of 
the function j{ 1 ' ' '(X; r,y) for each r G Z d , as A — > 0+. By the definition of Ji(A; r, y) after 
applying the Leibniz formula to its ([rf/2] — l)-th derivative we obtain 



[d/a]-i / 1 N ( fc ) 



l"/ -J ^ / | \ {K) 

jr^\X;r,y)= g ^M {^ ) - G *(0,y) J ^* 



«^-i-*)( A; r). 



Analyzing every summand in the last sum with the help of relations ( l59]) -( l6"Tj) and (172]) . ( l73|) 
we conclude that the first summand gives the main contribution to the asymptotic behavior of 
Jf l (A; x, y), A — > 0+. Indeed, the first summand has the order of growth 1/VX for uneven 

d and has the order 1/A for even d, respectively. On the other hand, the rest summands are 
0(1) for uneven d and are 0(ln 2 A) for even d, respectively. Therefore, for each r G Z d , as 
A — > 0+, one gets 

(-l^W^W.y)- { (gg^ gg^ ' if d isuneven ' (76) 

I (d-2)(Gg(o,o)-Gg(o,j,))A' lf rf is even. 

To estimate the asymptotic behavior of the functions J2 (X;r,y) and J3 (X;r,y), 

as A — )• 0+, we use the Leibniz formula once again, namely, 

4^-)(A; r, ») = '|f Cfo_, ( ffl>(0 °ffg (0|g) ) " (^(A; y - r) - *,(* „))«*-'>, 

I* 2 '- 1 / /-■ ni ...1 \ (*) 



4^-»(A; r , ,) = £ C &/2M ( ffl(0 °; ( _°'g ( o, g ) ) <*«* °> - *'<* '»' 



^([d/2]-l-fc) 

v ''--- J \cXia.o)--CXu). ll )J —•»< — < 

fc=0 

By virtue of (I5"9l) - (I6T1) . (I74|) and (I7S1) at the right-hand sides of the last equalities there are only 
the first summands which grow as o(l/vA) for uneven d and as o(l/A) for even d, respectively, 
as A — > 0+. The rest summands are 0(1) for uneven d and are o(ln A) for even d. Thus, on 
account of (1761) for each r G Z d , as A — > 0+, we come to the desired estimates 



4 d/2] - l \^y) = o (j?^^(A;0,y)) , J^^y) = o (j?^^(A;0,y)) . (77) 

Let us study the asymptotic growth of the function j] (A; y) as A — > 0+. The definition 

of J^A; ?/) and the Leibniz formula entail 

[<*/2]-l / , \(k) 



Due to relations (jo9l) - (lbT|) . (1721) and (1731) we note that the first summand at the last sum is 
o(l/vA) for uneven d and is o(l/A) for even d, respectively, as A — > 0+. Meanwhile, the rest 
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summands have the order of growth 0(1) for uneven d and o(ln A) for even d, respectively. 
Hence, taking into account f J76|) we conclude that for each r e Z d 

jr 2 ^(A;,) = o(jr 2M) (A;r,y)), A^0 + . (78) 

We turn to estimating of the asymptotic behavior of the function Jg 1 ' '~ (A; y), as A — > 0+. 
The definition of Js(A; y) and the Leibniz formula imply that 

[<*/2]-l X . X (fe) 

^""<* »> = E c ^- ( G>( o,o)- G> (o,») ) «* 0) - *<<* ^^ 



k=0 



On account of relations (l59l)-( )6T]) and ( |63i) we see that the first summand at the last sum is 
o(l/vA) for uneven d and is o(l/A) for even d, respectively. All the other summands have the 
order of growth 0(1) for uneven d and o(lnA) for even d, respectively, as A — > 0+. By (!76|) it 
follows that for each r e Z d 

#/ 2 ^(A;,) = o(jr 2M) (A;r,,)), A^0 + . (79) 

The last step is the revealing the asymptotic properties, as A — > 0+, of the function 
Jq (X;r,y). By the definition of Je(A;r, y) together with the Leibniz formula one can 

write 

T ([d/2]-i) (x ,_ y- rk ( 1 \ / OA(0,r) \ 

Recalling relations (155J1 . ( 157)) . ( 172")) and ( 175)) . we deduce that if g? is uneven then all the summands 
at the last sum are 0(1), as A — > 0+. However, if d is even then the first and the last summands 
have the order of growth 0(ln A) whereas the others have the order of growth 0(1), as A — > 0+. 
In any case in view of (176)) the following formula holds true 



4 m - l) (A; r,y) = o (j[ W ^ (A; r, „)) , A -► + . (80) 

Combining results (J551) - (J55|) and (ITS)) - ()8"0l we conclude that 

f°° -XtM/2]-l ( rr (^ rr mW - 2 7d (p d (y - x) + gg(g) - p d (y)) 



X 



>/7r/vA, if d is uneven, 
1/A, if d is even, 



/*oo /*oo 

/ e" A< tt^- 1 (^o, s/ ,o(oo) - H 0tyfi (t)) dt ~ / e- xt m-\H y!y , (oc) - H y , y>0 {t)) dt 
Jo Jo 

27d J y/n/^/X, if <J is uneven, 



i(d-2)(G (0,0) + G (0,!/)) 2 I 1/A, if d is even 



«) 



as A — > 0+. If Cd(x,y) > for all x,y G Z d \{0}, s 7^ y (it is easily seen that Cd(0,y) and 
Orf(y, y) are always positive for d > 3) then application of Corollary 43 in [26J to the last 
asymptotic relations implies (ITTj) . 
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Thus, to complete the proof of Theorem [T] for d > 3 we only have to verify the positivity 
of Cd{x,y) for x, y G Z d \{0}, x ^ y. The main idea of the demonstration is the same as 
that exploited for proving the positivity of C^-, •) in Subsection 4.2. Namely, we assume the 
contrary, that is, for some d G N, x G Z d \{0,y} and y G Z d \{0} the relation Cd(x,y) = is 
valid. This is equivalent to the following relations 

G (0, 0) + G (0, y) = G (x, y) + G (0, x) & G (0, 0) - G (0, x) = G (x, y) - G (0, y). (83) 

Note that the formula for H y ^ x ^(oo) appearing in dSJ) and the definition of function Pd(-) allow 
us to write 

_ Gq(0, 0)(Gq(x, y) - Gq(0, y)) + G (0, y)(G (0, 0) - G (0, x)) 
y ' X ' 0[OO) Gl(0,0)-Gl(0,x) 

Combining (I8"3"j) and the last equality we see that 

(G (0, 0) - G (0, x))(G (0, 0) + G„(0, y)) 



H VtX ,o(oo) 



G§(0,0)-Gg(0,z) 



G (0, 0) + G (0, y) G (x, y) + G (0, x) 



G (0,0) + G (0,a;) G (0, 0) + G (0,x) 

However, by formula fl30|) one has H y>x>0 (oo) < H yjX (oo) and due to fl20l) it follows that 

G (x,y) + Go(0,x) G (x,y) 

G„(0,0) + G (0,x) £ G^O^) ** Go( °' 0) £ °° {X ' V) - 

The obtained contradiction completes the proof of fill I) and the proof of Theorem [T] as well. 

5 Proof of Theorem [2 

Firstly note that in view of Corollary [2] it suffices to establish Theorem |2] for z = 0. Recall that 
for a simple random walk on Z the transition intensities a(r,r + 1) = a(r,r — 1) = a/2 and 
a(r,r + k) — for all r G Z and fc G Z such that |/c| > 1. By definition of the function (f)(9) 
and by that of the constant 7! one gets (f)(9) = a(cos9 — 1), 9 G [— 7T, 7r], and 7! = l/y2a/K. 
Hence, equality (J6]) and Lemma H] imply that pi(r) = |r|, r G Z\{0}. Observe that the formula 
for H X} y jQ (oo) in ([7j) is valid for a simple random walk on Z as well. Therefore, we infer that 
H x ,y t o(oo) = for y < < x and x < < y, H x>y;0 (oo) = x/y for < x < y and y < x < 0, 
H X) y >0 (oo) = 1 for < j/ < a; and x < j/ < 0, H 0)y>0 (oo) = l/(2\y\) and H yi y j0 (oo) = l-l/(2\y\). 
Thus, for y < < x and x < < y one has H xy $(-) = 0, that is, relation (|TS|) is proved. 
For < y < x and x < y < due to (!30|) we conclude that H X! y }0 (oo) = 1 = H Xjy (oo) and 
H Xj o >y (oo) = 0. This is equivalent to relations H Xj o, y (-) = and H X}Vj0 (-) = H X}V (-). Since the 
asymptotic behavior of the function H Xjy (oo) — H x>y (t) is established in Lemma HJ relation ( TT2|) 
is also proved. Moreover, formula ( T4"3"j) is true for a simple random walk on Z and, consequently, 
applying Tauberian theorem (Theorem 4 in [T4], Ch.13, Sec. 5) we find the asymptotic behavior 
of H y ^ >0 (oo) — Hy tyt0 (t), &s t — > oo , with Ci(y,y) = l/(2y/an) > 0. Thus, relation f fl3|) is 
proved. So we have to investigate only the asymptotic properties of H Xiy)0 (oo) — H Xjy>0 (t) when 
x G [0, y), x G Z, that is, to verify relations ( !i~4l) and ( TTBT) . 
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To this end recall a well-known result (the gambler ruin problem) for embedded chain 
{S n , n G Z+}. Namely, P(0 < S k < y, < k < n\S(0) = x) < (1 - e ) n for some e G (0, 1) and 
x G [0,|/], i6Z (see, e.g., [21], Ch.l, Sec. 9). Bearing on this result we derive that 

P (t < r y , < oo, N(t) = n\ S(0) = x) 

n=0 

P (0 < S k < y, < k < n, N(t) = n\ S{0) = x) 

71=0 

Eoo 
P (0 < S k < y, < k < n\ S(0) = x) P (N(t) = n\ S(0) = x) 
n=0 

< y 00 (i- eo) n(ggy! = e ^. 

It follows that relations (TT4|) and ( 1T5|) are valid with e e (0, Eq). Theorem[2]is proved completely. 

6 Proof of Theorem [3 

On account of Corollary [2] it is sufficient to prove Theorem [3] for z = 0. Since T y> o = t~ + r 
a.s., one has H Xt y t0 (t) = H~ y0 * G(t), t > 0. This immediately implies H x ^ (oo) = H~ y0 (oo), 
x G Z d , y G Z d \{0}. Furthermore, by formula (1171) we get that for A > 

(tf-^oM^ H^, )(A) = (H x>y ,o^o)^ H x , y>0 )(\) - a-'H^oiX). (84) 

In view of ([9]) - ( TT3l) and Tauberian theorem (Theorem 4 in [T4], Ch.13, Sec. 5) it is easy to 
see that the first summand at the right-hand side of (1841) makes the main contribution to the 
asymptotic behavior of the left-hand side of (|84"|) . More exactly, this is true if d < 3, x G Z d , 
y G Z d \{0}, except for a simple random walk on Z in the cases 0<x<y, y<x<0, 

i<0<j/ort/<0<i. Therefore, (i^ (oo) -#^ )0 )(A) ~ (^^0(00) - H x>yj0 )(X), as 
A — ¥ 0+, and due to (J9"|)- (fl3"|) and Tauberian theorem (Theorem 4 in [14], Ch.13, Sec. 5) we 
prove Theorem [3] for the mentioned d, x and y. 

When the random walk on Z is simple and 0<i<t/or|/<x<0, the estimate of 
H~ y0 (oo) — H~ y0 (t), as t — > 00, is a direct consequence of Theorem [2] and the inequality 
H x ,y,o{oo) — H Xt y t0 (t) > H~ (oo) — H~ y0 (t) used in the proof of Theorem [TJ for d = 1. 
Moreover, for a simple random walk on Z when x < < y or y < < x we may assert that 
H~ y0 {t) = by virtue of the established identity H~ (oo) = H Xt y yQ (oc) and relation ffT6|) . 
Thus, for these cases Theorem [3] is also proved. 

To find the asymptotic behavior of H~ {po)— H~ (t), as t — > 00, for <i > 4, we differentiate 
[rf/2] — 1 times the both sides of equality (1541) w.r.t. A > 0. The asymptotic relation for the 
([rf/2] — l)-th derivative of the first summand in (18"4"|) is given by formulae (jSTj) and (I8"2"j) . Let 
us analyze the asymptotic properties of the ([d/2] — l)-th derivative of the second summand in 
Integration by parts gives 

/ e-WM- 1 dH Xiy , (t) = X e- x HWy- l (H Xty>0 (oo)-H Xty>0 (t))dt 

Jo Jo 

- l) r e- x HW^- 2 (H Xiyt0 (oo) - H x>y>0 (t)) dt. (85) 
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Obviously, the first summand in (|55|) is o (/ °° e~ Xk i^ d ^~ x {H x>y j Q {oo) - H X) y )0 (t)) dt), as A -» 0+. 
By formulae (ITTj) . flHTj) . (182]) as well as Tauberian theorem (Theorem 2 in JT3], Ch.13, Sec. 5) we 
infer that the second summand in ( I8"5"j) is also o (J e _At t' d / 2 l~ 1 (if :r)yj o( 00 ) — Hx,y,o(t)) dt) , as 
A —?• 0+. Thus, the ([d/2] — l)-th derivative of the second summand in ( I84p does not make a con- 
tribution to the asymptotic behavior of the ([d/2] — l)-th derivative of the left-hand side of (I84j) . 
Therefore, /fe-^H^o^) - H~ y>0 (t)) dt ~ J °° e^HW^\H x , y , (oo) - H x , y , (t)) dt 
where the asymptotic relation for the last expression is provided by formulae 1J8T1) and ( J82l) . 
Consequently, employing Corollary 43 in [26J we complete the proof of Theorem [3] when d > 4. 
Theorem [3] is proved. 

The author is grateful to Associate Professor E.B. Yarovaya for permanent attention and to 
Professor V.A. Vatutin for valuable remarks. Special thanks are to Professors I.Kourkova and 
G.Pages for invitation to LPMA UPMC. 
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